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Abstract 

In this paper we give all the details of the calculation that we presented in our 
previous paper arXiv:0908.0387 where the infrared structure of the MHV gluon 
amplitudes in the planar limit for = 4 super Yang-Mills theory was considered 
in the next-to-leading order of perturbation theory. Explicit cancellation of the 
infrared divergencies in properly defined inclusive cross-sections is demonstrated 
first in a toy model example of " conformal QED" and then in the real AA = 4 SYM 
theory. We give the full-length details both for the calculation of the real emission 
and for the diagrams with splitting in initial and final states. The finite parts for 
some inclusive differential cross-sections are presented in an analytical form. In 
general, contrary to the virtual corrections, they do not reveal any simple structure. 
An example of the finite part containing just the log functions is presented. The 
dependence of inclusive cross-section on the external scale related to the definition 
of asymptotic states is discussed. 

Keywords: Super Yang-Mills Theory, Infrared safe observables, Maximally helicity 
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1 Introduction 



In recent years remarkable progress in understanding the structure of the planar limito 
of the A/" = 4 SYM (supersymmetric Yang-Mills) theory has been achieved. In the pla- 
nar limit this theory seems to be integrable at the quantum level and its possible solu- 
tion would be the first example of a solvable nontrivial four- dimensional Quantum Field 
Theory. The objects which were in the spotlight starting from the AdS/CFT (Anti de 
Sitter/Conformal Field Theory) correspondence [Jj were the local operators, namely, the 
spectrum of their anomalous dimensions. They were calculated on the one hand side 
from the field theory approach |2] and, on the other hand, as energy levels of a string 
in the classical background [31 S] revealing a remarkable coincidence. This coincidence 
being part of the general conjecture suggests the way towards solution of the model at 
the quantum level. 



1.1 Scattering amplitudes at weak coupling 

Other quantities of interest are the so-called MHV|§ scattering amplitudes. It was realized 
long ago that in the planar limit the pure non-Abelian gauge theories they do have a truly 
simple structure [5]. In papers [SI 13 [8], the powerful tool for calculating the loop expansion 
for these amplitudes was suggested which allows one to calculate the loop contributions 
to the amplitudes without calculating the usual Feynman diagrams the number of which 
grows exponentially with the growth of the order of perturbation theory. Even greater 
simplification occurs in the case of the J\f = 4 SYM theory where the loop expansion takes 
extremely simple form in comparison with a less supersymmetric case [9] . 

To see the hidden symmetries of the MHV amplitudes, it is useful to consider the 
color-ordered amplitude defined through the group structure decomposition 

^a-^oop) ^ gn-2 ^ Tr(T'^p(^)...T»M.))40(p^(^^, (1.1) 

Y lU/l y perm 

where An is the physical amplitude. An are the partial color-ordered amplitudes, T'^^^^ are 
the generators of the gauge group SU{Nc), api^i) is the color index of the p{i)-th. external 
particle, and is its momentum. 

To be more precise, it was found that these amplitudes revealed the iterative structure 
which was first established in two loops l^lOj and then confirmed at the three loop level by 
Bern, Dixon and Smirnov, who formulated the ansatz [TT] for the all-loop n-point MHV 

^Defined as g —J' 0; iVc —>■ oo; A = g^Nc fixed 

^MHV (maximally helicity violating) amplitudes are the amplitudes where all particles are treated as 
outgoing and the net helicity is equal to n — 4 where n is the number of particles. For gluon amplitudes 
MHV amplitudes are defined as the amplitudes in which all but two gluons have positive helicities. 
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where E^'' vanishes as e ^ 0, C(') are some finite constants, and M^\le) is the le- 
regulated one-loop n-point amphtude. 

It is not surprising that the IR divergent parts of the amphtudes factorize and expo- 
nentiate [12]. What is less obvious is that it is also true for the finite part 



where 'jcuspid) = Y.ii^^Ylcusp so-called cusp anomalous dimension [13] and Go{g) = 
^lii^y^o^ is the second function (dependent on the IR regularization) which defines 
the IR structure of the amplitude. 

According to the BDS ansatz, the finite part of the amplitude is defined by the cusp 
anomalous dimension and a function of kinematic parameters specified at one-loop. For 
a four gluon amplitude one has 



The cusp anomalous dimension is a function of the gauge coupling, for which four terms 
of the weak coupling expansion [2] and three terms of the strong coupling expansion [3l IH 
[13] are known. Integrability from the both sides of the AdS / CFT correspondence leads 
to the all-order integral equation [TH] solution to which, being expanded in the coupling, 
reproduces both series [T6] . 

For n = 4, 5 the BDS ansatz goes through all checks, namely, the amplitudes were 
calculated up to four loops for four gluons [2j for divergent terms (see also [17] for checking 
at order 1/e) up to two loops for five gluons [18] and up to three loops in [19]. However, 
starting from n = 6 it fails. The first indication of the problem was the strong coupling 
calculation in the limit n ^ oo [20j where the authors compute the value of the amplitude 
for a particular kinematic configuration for a large number of gluons and find that the 
result disagrees with the exact value of the amplitude from the BDS formula. The second 
indication came also from this duality, namely, from the comparison of hexagonal light- 
like Wilson loop and finite part of the BDS ansatz for the six-gluon amplitude. It was 
found that the two expressions differ by a nontrivial function of the three (dual) confor- 
mally invariant variables [21]. The third indication appeared in [22] where the analytical 
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structure of the BDS ansatz was analyzed and starting from n = Q the Regge hmit factor- 
ization of the amphtude in some physical regions failed. Finally, it was shown by explicit 
two-loop calculation [23] that the BDS ansatz is not true and it needs to be modified by 
some unknown finite function, which is an open and intriguing problem. However, from 
the two-loop calculation for the six-point amplitude and hexagonal light-like Wilson 
loop [21] it was shown that the gluon amplitude/ Wilson loop duality [25] is still valid. 

1.2 Strong coupling dual of amplitudes, light-like Wilson loops 
and dual conformal invariance 

In [25], the authors defined the prescription for calculating the amplitudes at strong 
coupling. It happens that in leading order the amplitude is given by the light-like Wilson 
loop living on the boundary of dual AdS space 



where S^i denotes the classical action of classical solution of the string worldsheet equa- 
tions in Euclidean space-time, which is proportional to the area of the string world-sheet. 

After this in [25j it was conjectured that duality between light-like Wilson loops and 
MHV scattering amplitudes is valid at any coupling, which was proved for n-point MHV 
amplitudes at one loop [27J and forn = 6 at two loops [24J (for more details and references 
see the review [2.8J). 

Due to the cusps the light- like Wilson loop is UV divergent; however, this divergency 
is under control, namely one can write the divergent factor in all orders in the coupling 
governed by two functions, one of them being the cusp anomalous dimension mentioned 
above. This allows one to define the finite parts for both the Wilson loop with n cusps 
and the n-point MHV amplitude which, according to DKS conjecture [25j, are equal to 
each other 



In [2n] the notion of dual superconformal symmetry was introduced, which is con- 
formal invariance acting in momentum space. What is important, this symmetry has 
a non-Lagrangian nature. After this in [301 EI] the fermionic T-duality was suggested 
which maps the dual superconformal symmetry of the original theory to the ordinary 
superconformal symmetry of the dual model. 

For a Wilson loop the conformal invariance is broken due to the cusps, but one can 
write the anomalous Ward identities which allows one to find the finite parts of the Wilson 
loop with n = 4 and n = 5 cusps exactly [32] 




Mn ~ exp[-S'|] = exp[— (Area)ci], 
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where the connection between the momentum space and its dual = — xf^^ = 
is used. This equation uniquely fixes the finite parts of the Wilson loop with n = 4 and 
n = 5 cusps; however, starting from n = Q more input is needed since the finite part of 
the Wilson loop in this case can be a function of the three conformal invariant variables. 
Hopefully, one can find hidden symmetries which fix the finite part for any n [33]. 

It is not clear how to derive this duality from the field theory point of view, and also 
how to extend it to the NMHV cas^. At one-loop one can show that finite part of the 
so-called two mass easy box which governs the finite function of MHV amplitudes could 
be directly mapped to Wilson loop diagrams through a simple change of variables in the 
space of Feynman parameters and also through the connection between scalar integrals 
in different dimensions 1341. 



1.3 Infrared-safe observables 

While all the UV divergences in A/" = 4 SYM are absent in scattering amplitudes the 
IR ones remain and are supposed to be canceled in properly defined quantities. By 
themselves the divergent amplitudes have no sense. Regularized expressions act like some 
kind of scaffolding which has to be removed to obtain eventual physical observables. It 
is these quantities that are the aim of our calculation. And though the Kinoshita-Lee- 
Nauenberg [35] theorem in principle tells us how to construct such quantities, explicit 
realization of this procedure is not simple and one can think of various possibilities. The 
well known example is a successful application to observables in QED |3S]. The other 
suggestion is to consider the so-called energy fiow functions defined in terms of the energy- 
momentum tensor correlators introduced earlier (see for example [37]) and considered in 
the weak coupling regime in pHl |39] and recently in the strong coupling regime in [IQ] . 
From our side we concentrated on inclusive cross-sections in hope that they reveal some 
factorization properties discovered in the regularized amplitudes. Similar questions were 
discussed in [4T], where the inclusive cross-sections like the IR safe observables based on 
on-shell formfactors in J\f=A SYM were constructed. 

To perform the procedure of cancellation of the IR divergences, one should have in 
mind that in conformal theory all the masses are zero and one has additional coUinear 
divergences which need special care. In this work we employ the method developed in the 
QCD parton model [121I1SI Sll 113 US] ■ It includes two main ingredients in the cancellation 
of infrared divergencies coming from the loops: emission of additional soft real quanta 
and redefinition of the asymptotic states resulting in the splitting terms governed by 
the kernels of the DGLAP equations [13 HH]. The latter ones take care of the coUinear 
divergences. 

Typical observables in QCD parton model calculations are inclusive jet cross-sections, 
where the total energy of scattered partons is not fixed since they are considered to be parts 

•^NMHV (next to maximally helicity violating) amplitudes are the amplitudes where all particles are 
treated as outgoing and the net helicity is equal to n — 6 where n is the number of particles. 
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of the scattered hadrons. In |13], the algorithm for extracting divergences was developed 
which allows one to cancel divergences and apply numerical methods for calculation of 
the finite part. In our paper we choose as our observables the inclusive cross-sections with 
fixed initial energy and get an analytical expression for the finite part of the differential 
cross-section. We do not assume any confinement and consider the scattering of the single 
parton based "coherent" states^ being the asymptotic states of conformal field theory. 

There are some attempts to deal with the divergences for the amplitudes themselves. 
For example in [35] a deformation of the free super conformal representation by contri- 
butions which change the number of external legs was proposed which looks similar to 
the procedure that we apply below considering the inclusive cross-sections. Acting along 
the same lines the authors of |50] have observed that the holomorphic anomaly [51] gives 
an extra modification of superconformal algebra for the tree level scattering amplitudes. 
They argued that superconformal symmetry survives regularization and introduced a new 
holomorphic anomaly friendly regularization to deal with the IR divergences. 

The paper is organized as follows. In Section [21 we consider general issues concerning 
the construction of the infrared-finite observables in the massless QFT. We discuss the IR 
and coUinear divergencies for the scattering amplitudes and the ways of their cancellation 
based on the Kinoshita-Lee-Nauenberg theorem. We introduce the notion of the measure- 
ment functions and discuss their properties. Then the concept of the splitting functions 
and splitting counterterms is outlined. We define the IR finite inclusive cross-sections 
which are the subject of calculations in the subsequent sections. 

Section [3] is devoted to the demonstration of the techniques discussed above in practice. 
In a toy model of "conformal QED" we consider the correction to the massless electron- 
quark scattering. We show how the IR and collinear divergences cancel and calculate 
analytically the remaining finite part of the differential cross-section. Due to absence of 
the identical particles in the final state this example turns out to be much simpler than 
gluon scattering in = 4 SYM and serves as a good warm-up exercise before going to 
U = A SYM. 

In Section IH we calculate the leading order PT (Perturbation Theory) correction to 
the gluon-gluon scattering inclusive cross-section. It includes the one-loop contribution 
to the 2 — )■ 2 scattering differential cross-section, the tree level 2 — )■ 3 scattering with 
the integration over the phase space of the fifth gluon and an account of the splitting of 
the initial and final states. We consider also the amplitudes with creation of pairs of the 
matter fields from the A/" = 4 supermultiplet. 

In Section [5l using the results of Section |4] we present the infrared finite results for the 
differential inclusive cross-section in A^ = 4 SYM theory for different physical setups. 

Section |6] contains discussion and concluding remarks. 

In appendices we present the technical details of our calculations. 

"^the squared perturbative amplitudes used in our calculation are summed over colors, so in this sense 
they are colorless and there are no contradiction with statements that cancellation of IR divergences 
occurs only for colorless objects. 
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2 Construction of Infra- Red Safe Observables 



The tree level matrix elements are finite and well defined in perturbation theory. Diver- 
gences appear when integrating over virtual loops or over phase space of real particles. So 
the first step is to choose a proper quantity which is finite in the lowest order of PT prior 
to calculation of radiative corrections. For example, the total elastic 2x2 cross-section 
is divergent, but differential cross-section is well defined. The choice of a proper quantity 
is performed by imposing conditions on the phase space. This can be achieved by intro- 
ducing the concept of the measurement function Sn, where n is the number of particles 
in the final state. It defines which physical quantity we are measuring. Typical examples 
are: a total cross-section, a differential cross-section, an n-jet cross-section, etc. In the 
case of 2 X n scattering the differential cross-section is given by 



da2^n 1 



- J \M2+n\'^d(j)n<Sn, 



where Sn is the measurement function and the n-particle phase space d(j)n is given by 

d<Pn = n^''b^)7Tr^(2vrA''(Pl+P2-p3-...-Pn+2). (2.8) 

fc=3 y'^'"') 

Here J is the flux factor, pi,p2 are the momenta of the incoming particles, p3,...,Pn+2 
are the momenta of the outgoing ones, M2+n is the matrix element of the corresponding 
process and we use the dimensional regularization with D = 4 — 2e. 
Then, for example, choosing the measurement function to be 



S2 = 6 {^Det - ^ 



13j 



one singles out the standard differential cross-section for the scattering of a third particle 
on a certain solid angle for the 2 — 2 process 



da2^ 



dn 



13 



M^\'d(f)2 S2. (2.9) 



If one wants to construct the IR finite quantity then, according to the Kinoshita-Lee- 
Nauenberg [35] theorem, it is not sufficient to consider the process with the fixed number 
of final particles. One has to include the processes of the same order of the perturbation 
theory with emission of extra soft quanta and integrate over their momenta. This leads 
to the notion of inclusive cross-section when one fixes some particles and integrate over 
all the others allowed by the conservation laws. 

When the number of particles increases, one has to specify the measurable quantity 
in a more accurate way and to distinguish the particle(s) in the final state. Thus, one 
can introduce the energy and angular resolution for the detector and cut the phase space 
so that the soft quanta with total energy below the threshold as well as all the particles 



8 



within the given sohd angle are included. This procedure requires the corresponding 
measurement functions and works well in QED but introduces explicit dependence on the 
energy and angular cutoff, thus violating conformal invariance. 

We adopt here different attitude without introducing any cutoffs but rather considering 
the inclusive cross-section with the emission of particles with all possible momenta allowed 
by kinematics. Having identical particles in the final state one has to specify which 
particles are detected by introducing some measurement function. For instance, one can 
detect the given particle scattering on a given angle while integrating over the phase space 
of the other particles. As it will be clear later in this case due to collinear divergences 
one still cannot avoid introducing some scale related to the definition of the asymptotic 
states of a theory. Below we show how it works in particular examples. 

To have the cancellation of all the IR divergences, according to the analysis of Ellis, 
Kunszt and Soper [13], the measurement functions for the processes with a different 
number of external particles have to obey the following conditions: 

Sn+i{...,Xp,...) =Sn{ ), A^O, (2.10) 

which reflects the insensitivity to the soft quanta, and 

Ap, (1 - A)p, ...) = ...) (2.11) 

here < A < 1. This condition expresses our insensibility to collinear quanta. 

It should be pointed out that in case of identical particles one has an additional 
problem when calculating the differential cross-sections: one has to specify the scattering 
angles and to choose the detectable particle. This requirement imposes further conditions 
on the phase space as will be shown below when considering the gluon scattering. 

The additional divergences appearing in the massless case which come from the in- 
tegration over angles rather than the modulus of momentum, as in the case of the IR 
divergences, are related to the collinearity of momenta of two particles. For this reason 
they are called the collinear divergences. To get the cancellation of all divergences, the 
observed cross-section should include besides the main process and emission of the soft 
quanta the process of emission of collinear particles with kinematically allowed absolute 
values of momenta. As we will see below, the leading IR divergences coming from the 
cross-section of the processes with the virtual loop correction and from the real emission 
of the soft quanta cancel. However, the total cancellation of divergencies does not happen. 
The remaining divergences in the form of a single pole have the collinear nature. For the 
cancellation of the remaining pole one has to properly define the initial (and final) states. 
The reason is that a massless particle can emit a collinear one which carries part of the 
initial momentum and in this case, it is impossible to distinguish one particle propagating 
with the speed of light from the two flying parallel. This is the common problem for any 
theory containing the interacting massless particles. 
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To deal with this problem, let us consider a particle in the initial state and introduce 
the notion of distribution of a particle with respect to the fraction of the carried momen- 
tum z: q{z). Then the zero-order distribution corresponds to q{z) = 6{1 — z) and the 
emission of a collinear particle leads to the splitting: the particle i carries the fraction 
of momentum equal to z, while the collinear particle j - (1 — 2;). The probability of this 
event is given by the so-called splitting functions Pij{z) [48]. In case of a particle in a 
final state, this corresponds to the fragmentation into a pair of particles i and j. In the 
lowest order of perturbation theory the distribution can be written in the form 

= *(i-^) + ^i(|J)'i:p«w. (2.12) 

where the scale Q'j, sometimes called factorization scale, defines the measure of collinearity 
of the emitted particles, i.e., it refers to the definition of the initial state. In fact, in the 
massless case one cannot define the initial state that contains just one particle, it exists 
together with the set of collinear particles forming a coherent state. 
This leads to the additional terms in the cross-section 

^] fdzY,P,,{z) daf°^-{zpuP2,P3,P^) + {Pi^P2), (2.13) 

referred hereafter as the initial splitting contributions or collinear counterterms. 

The same is true for the final states. The corresponding final state collinear countert- 
erms are ^ 

^) l'dxJ2P.,{x) rfaf ™(pi,p2,P3,P4) + (P3 ^P4). (2.14) 

Summarizing all the contributions we come to the following set of IR safe observables 
that we consider here 

d^Tbf = jdzi qi{zi, Jdz2 q2{z2, ^) H J ^Zi qi{zi, ^) x (2.15) 

xda^^'^{z,p,,Z2P2,...)Sn{{z})=g'NtY. (^f^j ^^f ™*^(5, t, Qj), 

where Pi-,P2 are the momenta of the initial particles, Pi are the momenta of the final 
particles, 5„ are the measurement functions which define the measurable quantity, qi are 
the initial and final state distributions. 

The above expression looks like the parton model cross-section. The difference is 
that in the parton model one uses both the parton distributions inside hadrons and 
fragmentation functions for final-state hadrons while here it belongs to the definition of 
the asymptotic states. 
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3 Toy model: "Conformal QED 



To illustrate the main ideas of the previous chapter, we study first a toy model example. 
Let us consider the electron-quark scattering and put all the masses equal to zero. We 
will be interested in the radiative corrections in the first order with respect to the strong 
coupling Us . The corresponding diagrams are shown in FigHJ We define the initial electron 




a) b) c) d) e) f) 



Figure 1: The process of electron-quark scattering in the first order in a^: ) the Born 
diagram, b)-d) the corrections due to the virtual gluons, e)-f) the corrections due to the 
real gluons 



and quark momenta as Pe = pi and Pg = p2 and the final ones as pe' = ps and Pq/ = pi, 
respectively. 

In the chosen process the UV divergences cancel due to the Ward identities so we are 
left only with the IR ones. To handle them, we use dimensional regularization. This 
situation exactly imitates the four- dimensional CFT's like the A/" = 4 SYM theory. 

Define the measurement function in the following way: 

where S±^h3 means that we detect the third particle with any helicity, i.e. we are interested 
in unpolarized differential cross-section. Here dilis = (i0i3(icos(6'i3)lj, 6*13 is the scattering 
angle of the particle with momentum pg with respect to the particle with momentum p^^ 
in the center of mass frame. In the leading order (LO) we have the well known text-book 
formula [52] 




(3.16) 



Born 

where E is the total energy of initial particles in the center of mass frame, and s,t, m are 
the standard Mandelstam variables. In the cm. frame s = E'^,t = —E'^/2{1 — c),u = 
-EV2(1 + c), C = COS^i3. 



^if to be more accurate in dimensional regularization we have dfi^^ ^ 
^'l3sm((/.l3)-2^rfcos(0l3)sm(0l3)-'^ D = A-2e. 
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The one loop correction coming from the diagrams with virtual gluon, Fig|T]b)-d), has 
the form 



'da2- 



da' 



2-s>2 



a 



13 / I— loop 



~2CF-^r^ (-T + - + 



13 / Bo 



Art 



-t 



(3.17) 



In order to avoid the transcendental numbers, we used the helpful definition of the angu- 
lar measure in the space of 4 — 2e dimensions and multiplied the standard expression by 
r(l — e)/(47^)^ Then the constants like •je, iog^An) and ({2) disappear from the inter- 
mediate expressions. Due to the cancellation of divergences in the final expressions, this 
redefinition does not influence the answer. 

Now, following the general prescription, we have to calculate the diagrams with emis- 
sion of real gluons, Fig|T]e) - f). For the measurement function of these processes we take 
<S3 = S2, then all requirements on S2 and 1S3 are satisfied trivially. Besides the squares 
of each of the diagrams one should also take into account the interference term. After 
contracting all the indices the phase integral takes the form (we denote the quark mo- 
mentum as p4 — and the gluon momentum as k and keep the standard notation for the 



Mandelstam variables s = {pi + P2) , t = {pi — p^) 



u 



(P2 - Ps) 



' da2. 



^ Jd^PsS^pl) j^S^{e)6^{{p.-kf)Sm 



dn 



13 



Born 

|M|2 = 

Mo 
Ml 



(27r 



P4=P1+P2-P3' 



4 t{s + t + u) 
4s — 8pik — 4p2k + 



(3.18) 



-8{pikf + 4(2s + t)pik - (3s2 + t^ + u^ + 2st) 



:^(s + u) + 8piA;+8p2fc-f 
■-A{s + t + u)- Ap2k - 



P2k 

%{pikf -4:{s+t+u)pik+2{s+t+uf -2{u + s)t 



[s + t + u] 



P2k 



P2k 

[s + t + u- 2p2kf 
P2k 



It is useful to pass to the spherical coordinates and use the cm. frame, 
integration over the phase volume the result can be represented as 



After the 



'da' 



2-i.3 



' da^ 



2^2 



dVt 



13 / Born 



dn 



13 / Bo 



2C, 



+ C.^^i^] (^) i^ + h) + Oie 



^2 An 



s 



An 



2 3 



t 



e 



(3.19) 



where the functions /i and /2 in the cm. frame are 

(l-c)(c3 + 5c2-3c+5)log(^)-(c-l)2(c+l)(c-ll)/4 



fi 



;i-c)2(i + c) 



(3.20) 



12 



2 



(1 - c){c^ + Sc^ - 3c + 5) log^(^— ^) 



+ ^(1 - c)(3c=^ + 15c2 + 77c-31) log(^^) + (l + cY{^ + 2c+h)u^ 
-12(9c2 + 2c+5)Lz2(^^) + ^(l - c)(l + c)(5c'-42c-23) 



(3.21) 



As one can see from comparison of the cross-sections of the processes with virtual 
fl3.17p and real gluons fl3.19p . in the sum the virtual part completely cancels and the 
second order pole disappears. However, the total cancellation of divergences does not 
happen. The remaining divergences in the form of a single pole have a coUinear nature. 
As was already mentioned, for their cancellation one has to define properly the initial 
states. 

Introducing the distribution function for initial quark state one gets the additional 
contribution j53| to the cross-section which looks likqj 



' do 2^2\ 1 «s 




d^P,,i^)\7^] (^(Pi,^P2) , (3.22) 



\^^13/5piit ^27r^u W// \ ""13 J Born 

where the Born cross-section is given by fl3.16p with the replacement of the initial quark 
momentum p2 by P2Z- This means that one should keep the total energy but replace the 
Mandelstam variables m according to eq. fl4.64p (see below). One should also keep the 
e term which gives contributions to the finite part. The splitting function Pqq{z) [IS] here 




(3.23) 

Notice also the change of momenta conservation condition which now looks like pi + zp2 — 
P'i~ Pa = 0. This gives an additional factor of Az'"^"^ / {1 + z — c(l — z)Y'^^~'^\ 

One might also have a contribution from the final state counterterm; however, since in 
this case, according to fl2.14p . the cross-section does not depend on the fraction z, one has 
to integrate only the splitting function Pqq{z). And this integral equals zero due to the 
requirement of conservation of the number of quarks. Therefore, one has no contribution 
from the final state splitting. It will not be the case for the gluon scattering cross-section 
considered below. 

The factorization scale Q'j is an arbitrary quantity associated with the quark distribu- 
tion function which may depend on z. It is quite natural to choose the factorization scale 
equal to the characteristic scale of the process of interest. Thus, in our case this choice 
corresponds to = —t, where t is the Mandelstam parameter t for the process where 



^We put here (^q^^ inside the integration over z since in general one may consider to be a functic 
of z. 
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P2 is replaced by p2Z. One has t = (see eq. (l4.64p below). Substituting this 

value of into fl3.22p leads to the following result: 



'da 



2^2 



13 



Split 



oP' as 



s 



(3.24) 



where 

h = 



^2\2 



2(1 -c)(c^ + c^ -33c + 7) log 



,1 -c. 



+ 12(9c2 + 2c + 5)Li2(^^) 



c r c- 



2c+ 5)7r2 - -(1 - c)(l + c)(llc2 - 19) 



2 

(3.25) 



Gathering all pieces together we finally obtain the IR finite answer in the NLO order 
of PT: 



da 



'da- 



2-^2 



dn 



13/ IR safe 



dn 



+ 



'da' 



2^2 



' da- 



2->3 



13 



a 



c2 + 2c + 5 



Bom 



as 



, dni3 

Cf 



l—loop 



dn 



+ 



'da 



2^2 



13 



Born 



2^2 [ (1 -c)2 27r(l -c)(l + c)^ 



[c^ + Sc^ - 3c + 5) log^ 



dni3 
,1 -c. 



Split 



-(7c^ + 19c^ - 55c - 3) log(— — ) - (1 + c)(3c^ + 21c + 2) 



(3.26) 



This expression is the final answer for the cross-section of the physical process of the 
electron- quark scattering where the initial and the final state include the soft and coUinear 
gluons. It includes also the definition of the initial state and can be recalculated for the 
alternative choice of the factorization scale similar to what happens to the ultraviolet scale 
which defines the coupling constant. Thus, we practically deal with the scattering not of 
individual particles but rather with coherent states with a fixed total momentum. This 
process contrary to the scattering of individual massless quanta has a physical meaning. 
The drawback is the dependence on Q/ which reflects the definition of the asymptotic 
state. This dependence explicitly violates the conformal invariance. 



4 Calculation of the Inclusive Cross-sections in J\f = A 
SYM theory 

Consider now the gluon scattering in the J\f = A SYM theory. Our aim is to evaluate the 
NLO correction to the inclusive differential polarized cross-section in the weak coupling 
limit in planar limit in analytical form and to trace the cancellation of the IR divergences. 

We start with the tree level 2 — )■ 2 MHV scattering amplitude with two incoming 
positively polarized gluons and two outgoing positively polarized gluons and consider the 



14 



differential cross-section da2^2{g^g^ g^9^)/dO. as a function of tlie scattering solid 

angle. Treating all the particles as outgoing this amplitude is denoted as ( MHV 

amplitude. At tree level the differential cross-section is given by 

1 y {tree) 1 2 



>2 



|A^rf52, (4.27) 



where J is a flux factor, in our case J = and the phase volume of the two-particle 
process (we use the FDH version of the dimensional reduction, see [5^ for details) is 

#2 = -^^^^pn (2^)i^-i (^^^ ^ (P1+P2-P3-P4), (4.28) 

and Sn {n = 2) in this particular case is 

S2 = 6+,h,S''-\^Det-^i3), (4.29) 

where 6^~'^{Qj:)et — ^13) means that our observable is the differential cross-section and 
5+^/13 indicates that we detect a particle with positive helicity. 

The squared matrix element is obtained from the color-ordered amplitudes via sum- 
mation 

|_^(*ree)|2 ^ ^4^2(^2 _ ^ | (pi , , P.(3) , P.(4) ) P , (4.30) 

where P„ is the set of the permutations of n objects (n = 3 in this case), so that in our 
case [HI [55] (see also appendix A for details) 

|^(i,ee)(--++)|2 ^ ^4^2^iv2-l) ^ , (4.31) 

CTgP3 'Sl<j{2) Sa{2)a{3) Sa{3)a{4) 5^(4)1 

where we use the notation Sij = {pi + pj)"^. The corresponding Feynman diagrams are 





Figure 2: Tree- level diagrams for the color-ordered MHV amplitudes 
shown in Figl2j 
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Within the dimensional regularization (reduction) the cross-section in the planar limit 
looks like 

V J {tree) 



Q?Nl I S" 



s 



(4.32) 



where s, t, u are the Mandelstam variables, E is the total energy in the center of mass 
frame, and a = g'^Nd^'K. So in the center of mass frame the cross-section can be rewritten 
as: 

''''^2-^2 



'da,^,V''^^^ a2iV2 4(3 + c2) //i^^ 



dV^^ 



(l-c 



2\2 



(4.33) 



/ {tree) 

where c = cos^is. The next step is to calculate the NLO corrections. 



4.1 Virtual part 








+ permutations 



I 



- fW 
I 



^A> — ' — *\A/ 



Figure 3: The one- loop diagrams for the color-ordered MHV amplitude in the A/" = 4 
SYM theory. Particles running inside the loop include all the members of the J\f = 4 
supermultiplet. The solid and dashed lines correspond to the fermion and scalar particles, 
respectively. 



To get the one-loop contribution to the differential cross- sect ion, one has to consider 
the diagrams shown in Figj3l We use the already known one loop contribution to the 
color-ordered amplitude [9] 



M, 



{l-loop) 



^{l-loop) I j^tree) 



s,t). 
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where ^°°^\s,t) is the scalar box diagram 



r(l-/oop) 



2 r(i-6)^ 1 

stV{l-2ey e2 



+ ■ 



TT 



0(e). 



The square of the matrix element summed over colors 



\M 



(i-;oop)|2^ ^ (4*^^^Ui'-'"°"^* + c.c.) 



colors 



has the form 



(l-ioop)( ++)|2 

4 I 



(4.34) 



s,t) -\ — t-^sulji' ^°°^\s,u] 



which gives the one-loop contribution to the cross-section in the planar limit 



t,u) 



-^2 



13 



+ (( 



(— ++) 

I 

virt 

r + 



2^2 



,2 



a 



47r s'^t'^u^ 
-u / 



4 ((^r+(— )V 



— ti 



16 

+—n\s^+t^+u^)+A{u^ log 



+ t^log^ 



' ^+^^log^(-)) 



(4.35) 



Rewriting this expression in the center of mass frame we have: 



'da- 



2->2 



2e 



a 



dn 



13 



virt 



5 - 2c + 
H — ^TT^ log I 



1 + c, 



16 3 + c2 



4 /5-|-2c + c2 



-,2^2 



(1 -C2)2 



16(3 + c2)7r2 

3(1-C2)2 " (1 



16 



+ - 



log( 



:i-c 



,2A2 



log I 



1 -c. 



c ^ 1 + c 



. (4.36) 



It should be stressed that due to the conformal invariance of the A/" = 4 SYM theory 
at the quantum level there are no UV divergences in fl4.36p and all divergences have the 
IR soft or collinear nature. They have to be canceled in properly defined observables. 
Note also the simplicity of the finite part which is a consequence of symmetries of A/" = 4 
SYM and the fact that all the terms in fl4.36p have the same transcendentality p3| 156] . 



4.2 Real emission 

The next step, as in the toy model considered above, is the calculation of the amplitude 
with three outgoing particles. Here we have to define the process we are interested in. 
There are several possibilities. 



17 



1. Three gluons with positive hehcities: g^g^ — ?■ g^g^g^ ■ This is the MHV amphtude; 

2. Two gluons with positive hehcities and the third one with negative hehcity: g^g^ — )■ 

This is the anti-MHV amphtude; 

3. One of three final particles is the gluon with positive helicity and the rest is the 
quark-antiquark paii@: g^g^ — t- g^q^q^ or g^g^ — t- g^q^q^ . This is an anti-MHV 
amplitude; 

4. One of three final particles is the gluon with positive helicity and the rest are two 
scalars: g^g^ — ?■ g^^h.. This is an anti-MHV amplitude. 

The corresponding diagrams are shown in FigJH 



Figure 4: The tree diagrams with three outgoing particles for the color-ordered ampli- 
tudes. Permutations are not shown. 

If one fixes one gluon with positive helicity scattered at angle 6 and sums over all the 
other particles then all the processes mentioned above contribute. In the case when one 
fixes two gluons with positive helicity and looks for the rest, only the first two options are 
allowed. 

^There is also a g^g^ — >> g^g^g^ helicity configuration. The partial amplitudes for both i cases 
where the additional gluon with negative helicity is the i second or the third gluon in the final state are 
equal. We will use i the ( h H — ) notation for both of them. 

®The J\f — 4: supermultiplet consists of a gluon g, 4 fermions ("quarks") and 6 real scalars A"^^; A 
and B are S'C/(4)fl indices, A is an antisymmetric tensor. It is implied that all squared amplitudes with 
quarks and scalars are summed over these indices. 
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The cross-section of these processes can be written as 
where ci^s is the three-particle phase volume 

?i +P2 -P3 -P4 -Ps), (4.38) 



d^'psS+jpl) d^p,5+ipl) d^p,5+ipl) ^ ^ 
(2vr)^-i (2vr)^-i {2n)^-^ ^^"^ ^ 



and is the measurement function which constraints the phase space and defines the 
particular observable. 

The squared matrix element is expressed through the amplitudes as before 

\Mt'\p„...,p,)\' = g'N!{N'^-l) Y: \At'\pi,P.w,...,P.(.))\'. (4.39) 

creP4 

For the processes mentioned above one has the following expressions for the matrix 
elements: [^: 

1. |7Wf = /iV3(iV2 - 1) ^ '-^ . (4.40) 

<^eP4 Si^il) ■Sa(l)a(2) Sa(2M3) S<T(3)a(4) Sa(4)l 

Since there are three identical gluons with positive helicity in the final state one has to 
define which ones are detected. In case of one detectable particle, one can choose the 
fastest one; in case of two, the two fastest ones. The measurement function for detecting 
only one gluon with momentum p^ with positive helicity can be written as 

5^-+++^'^ = hM^iPs > pIMpI > pDS^'-'innet - fii3); (4.41) 
and for detecting of two gluons with positive helicities as 

St'^^^^'' = 5^,,,5+,,MpI > pIMpI > pDS^'-'inoet - (4.42) 

where we detect the 3-rd and the 4-th gluons. Analogous measurement function would 
appear if we would like to detect the 3-rd and the 5-th gluons. 

2. \Mt""^^-^^-'^\^ = g^Nl{Nl - 1) E — ; (4-43) 

<^eP4 Sla{l) Sa{l)a{2) ■5<7{2)a(3) S<x{3)a(4) ^^(4)1 

The measurement function for detecting one gluon with positive helicity and momentum 
Ps is given by 

51--+++^'' = 5^,hMpl > pl))5^-\nnet - n,,); (4.44) 



^It is implied that all squared amplitudes with quarks and scalars are summed over SU{A)j^ indices. 
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and for detecting of two gluons with positive helicity by 

^(-++-),2 ^ s^^^^s^,^^^^^6^~2^noet - ^u). (4.45) 

3. \Mt""^^-^''^\' = g^'NliNl -l)Y. ^34535(4 + ^15) . (4 

aeP4 ^Ml) ■5^(l)'^(2) ■5<7(2)<7(3) Sa(3)<7(4) 5^(4)1 

The measurement function in this case is simple since we have only one gluon in the final 
state 

5^-+^^) = 5+,H.,6''-\nDet - ^13). (4.47) 

4. |_y^(*-e)(-+AA)|2 ^ ^6^3^^2 _ ^) ^ 44 _ ^4 43) 

ae Pi -^l-^Cl) ■^'^(1)^(2) ■^'7(2)(t(3) ■Sct(3)(t(4) 5^(4)1 

The measurement function is given by the same formula (14.471) as in the previous case. 

One can check that the measurement functions written above satisfy the IR and 
collinear limit conditions f l2.1U|2.1ip . Indeed, one has 

1. P5 0, |p4| = |p3|: 

^3(p3, p4, 0) ^ 0(j9° - vDd'^-^n - n,), 

2. p3 = -P,P4 = xP,p5 = (l-x)P: 

^3(p3,p4,p5) ^ e(i - x)Q{x)5''-\n - n,). 

The latter 6'-functions give < a; < 1 restricting the fraction of momenta in a natural 
way. 

Choosing the fastest momentum one has to have in mind the conservation of momen- 
tum and energy 

P3 + P4 + P5 = o, |p3| + Ip4| + IpsI = -E- 

This means that the three momenta form a triangle with the perimeter equal to E. Hence 
the requirement that, say, the third particle is the fastest one means that > E/3. 
Therefore, to simplify the integration, in what follows we choose the universal measure- 
ment function 

53(P3,P4,P5) = e(p° - ^E)S''-\Qnet - ^3), (4.49) 

where we take 6 = 1/3 in the case of identical particles and 5 = 1 in the other cases. Thus, 
the registration of one fastest gluon corresponds to 5 = 1/3 for the MHV and anti-MHV 
amplitudes and 6 = 1 for the matter-antimatter amplitude, while the registration of two 
fastest gluons corresponds to 5 = 1/3 for the MHV amplitude and 6 = 1 for the anti-MHV 
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amplitudj^. In what follows we keep the value of 5 arbitrarily and show that the IR and 
collinear divergences cancel in observables for any value of 6. We omit the details of the 
calculation, which can be found in Appendix B, and present here only the divergent parts 
of the calculated objects. All the finite parts can be found in Appendix D. 

With these definitions the contributions to the 2 — )■ 3 cross-sections from the ampli- 
tudes that are listed above are 

1. Real Emission (MHV) 



'da- 



2->3 



(—+++) 



a^N^ //i^y^a f 8 (3 + c^) 



dQ 



1 

+- 
e 



13 / Real 
2 



logi 



1-c. 2 



1 + 2 
-|-Finite part} ; 



-) + 



:i-c: 



J TT [ (1 — c^y 

, A + c, 165(25-3) 12(3+c2), A-5 
- log H ^ 1 ^ log' 



(4.50) 



;i-c2)2(l-5)2 



1-c 



6 



Notice the singularity in the limit 5 — )■ 1. 
2. Real Emission (anti-MHV) 

/da2-^V-^^-^ ^ /fj^V' a / 1 8(3 + c^) 1 



13 / Real 
+ 



12(c2 + 3)log5 



:i-c 



3\2 



64(12c2 + 17) 
3(1 -c2)3 



+ 



26 



'1-c 



.2\2 



-(5 + 3c^)6^ - (c^ + 19)5 + 2(5c2 + 43) 



+ 



'2(3c2 - 24C + 85) , ,1-c, 8(c2-6c + 21), l+5-(l-5)c, 
log(^^) - T, \n , N-^ t; ) 



(l-c)(l + c)3 2 
32(c2-4c + 7) 



(1 -c)(l + c)3 

32(2 - c) 



+ 



(l + c)3(l-c)(l + 5-c(l-5)) (l + c)3(l + 5-c(l-5))2 
64(1 - c) 



3(l + c)3(l + 5-c(l-5))= 



+ c 4-> -c 



+ Finite part > ; 



(4.51) 



Contrary to the MHV case the limit 5 — ?■ 1 is regular here and greatly simplifies the 
final result. 

3. Fermions (for 4 fermions in adjoint representation of SU{Ne)) 



' da- 



2^3 



( +9?) 



2(3 -c)^ 

:i-c){i + cy 



\ s 



1 / 1 ~ c. 



2(3 + c] 



16 

e 



;i-c)3(i + c) 



(79 + 25c^; 
3(1 -c2)2 



(4.52) 



+ Finite part 



I'^These are not precisely the needed requirements but are pretty close to them. Fulfillment of the 
exact requirements of the fastest particles is technically more involved but does not changed the picture. 
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4. Scalars (for 6 scalars in adjoint representation of SU{Nc)) 











'_8 







2(10 + 7c^ 



3(5 + c), ,1 + c. 




(4.53) 



In the last two expressions we chose the parameter 5=1 since there are no identical 
particles in these cases and there is no need to restrict the phase space. Note also the 
absence of the second order pole in e which means that there is no IR soft divergency here 
but only a collinear one. 



4.3 Splitting 

Now we have to deal with an additional 1/e pole coming from the collinear divergences. 
As one can see from the toy model example, taking into account emission of additional 
quanta in the initial and final states allows one to cancel the IR divergences (double 
poles in e) but leaves the single poles originating from collinear ones. Indeed, as it has 
been discussed earlier, the asymptotic states (both the initial and final ones) are not well 
defined since a massless quantum can split into two parallel ones indistinguishable from 
the original. To take this into account, we introduce the distribution of the initial and 
final particle (gluon or any other member of the = 4 SYM supermultiplet) with respect 
to the fraction of the carried momentum z: qi{z , Q'j / n"^) . Also, one has to keep in mind 
that the particles in this case are polarized. The corresponding Feynman diagrams are 
shown in FigJSl 

Additional contributions from collinear particles in the initial or final states to the 
inclusive gluon cross-section (the collinear counterterms) have the following form, respec- 
tively: 

dcTs^la™* = tJ-- ( ^ ) E / dzJ2 Pgi{z)da2^2izpi,pj,p:i,p4)S2^^'''''\z), (AM) 

Z'K6\LJj/ • 2; i^j l=g,q,A 



da^^J^i'"" = ( y-^] da2-,2{pi,P2,P3,P4) f dz J2 Pgi{z)Sf'^'''{z). (4.55) 



Having particles with different helicities we have the following set of collinear coun- 
terterms (We use here slightly different notation for the splitting functions indicating 
explicitly all three particles like Pfinijin2(^) avoid confusion.) 

1. Initial state splitting MHV amplitude ( 1- + +) 



da' 



(—+++) .,1 / ,,2\^ .1 /J„ \ (— ++) 



^^13 J InSplit 




^ ^ - iiM W o P.- ..W ^ fzp„p,,p,,p^st-;^-\z) 

(4.56) 
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e) f) 9) h) 



Figure 5: The initial and final particle splitting diagrams: a) the initial MHV amplitude, 
b)-c) the initial anti-MHV amplitudes, d)-f) the initial matter amplitudes, g)-h) the final 
MHV and anti-MHV amplitudes. Permutations are not shown. 



Final state splitting MHV amplitude ( 1" + +) 

(4.57) 

2. Initial state splitting anti-MHV amplitude ( h H — ) 

+ {Pi^P2), (4.58) 
Final state splitting anti-MHV amplitude ( h H — ) 

(tf - S)' (tf ) '----^r- ^»-<^)^r- -''^)- 

(4.59) 

One has also the coUinear counterterms containing the other members of the J\f — A 
supermultiplet 
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3. Initial state splitting into a fermion-antifermion pair 



^^13 ) insplit 27r e \Q)) Jo 



2 PfA^) (^\\Zlp2,P3,pdSt-:f\z) 



1'"^' ' \ dQ 



13 



+2 (^\\Zlp2,P3,P,)Str^f\ 



13 



+ (pi^P2), (4.60) 



4. Initial state splitting into a scalar pair 



'13 



+ (Pi^P2), (4.61) 

where Uf and is the number of fermions and scalars, respectively. One should put 
Uf = 4 and = 6 in our case. 

The explicit form of the Born cross-sections and the splitting functions Pjk{z) can be 
found in Appendices A and C, respectively. Notice that when changing the momentum 
Pi — )■ zpi one has to modify the Mandelstam variables according to eq. fl4.64p and take into 
account the additional factor from the phase space in full analogy with the QED case (see 
the comment after eq. fl3.22l) ). 

Note that there is no final state splitting counterterms for fermions and scalars. The 
reason is that one has to take into account only those final splittings where the original 
state (gluon in our case) survives with momentum multiplied by fraction z. 

The measurement functions here are the same as in the case of real emission but 
depend now on fraction z and restrict the integration region over z. They take the form 

Sf^\z) = (5+,,3(5^-2(fi - ni3)e{z - z™„), (4.62) 

for detecting of one gluon and 

-,h3 



for detecting of two gluons. 

The values of Zmin can be calculated from the requirement > {l — 6)E/2 in different 
kinematics. Indeed, for the initial splitting process one has to change the momentum of 
in-going particle, for example, pi to zpi which gives in the cm. frame 

S — zs, 

^ ~^ l + 2-c(l-z)^' 
^ ~^ l + 2-c(l-z)'"' 

Pi ^ TTI^f • (4-64) 
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At the same time, for the final splitting one has to substitute z^. This leads to the 
values of z^in, respectively, 



min 



{l-S)il-c) f,^ 



l + 5-c(l-5)' 



= (1-5). 



(4.65) 



Taking into account the splitting of the initial states and the fragmentation of the 
final states we get the following contribution to the inclusive cross-sections: 

1. The initial and final splitting for the MHV amplitude 
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+ Finite part > , 
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(4.66) 



(—+++) 

^^^13 JfuSpIU 



14(cM_3) 1-5 " 

QjJ [ "7(1-C2)2 



fi \ a 



(4.67) 



2. The initial and final splitting for the anti-MHV amplitude 
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(4.68) 
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■ (4.69) 



3. The initial sphtting for the fermion final states (5 = 1) 
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(4.70) 



+ Finite partj ; 
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4. The initial splitting for the scalar final states {6 = 1) 



'7 \ (— +AA) 




(4.71) 



5 IR Safe Observables in A/" = 4 SYM 

In the NLO there are two sets of amplitudes, namely, the MHV and anti-MHV amplitudes 
which contribute to the observables. The leading order 4-gluon amplitude is both MHV 
and anti-MHV and we split it into two parts. Then one can construct three types of 
infrared finite quantities in the NLO of perturbation theory, namely, 

• pure gluonic MHV amplitude 



^MHV _ 1 / da2~^\ ^ / da2^\ ^ ^^^^ / da2-^\ ^ ^^^^ f (i(T2_>2 V 



(5.72) 



pure gluonic anti-MHV amplitude 



2W^13/y,rt W^13 /iieaZ \ ^^^^ J InSpHt \ ^^^^ ) FnSplit 

(5.73) 

• anti-MHV amplitude with fermions or scalars forming the full A/" = 4 super mult iplet 



( + , qq+AA) / , s ( + , gg+AA) 

Real \ '^^13 / inSplit 



^Matter ^ / ' + ( ^ 1 • (5.74) 



We would like to stress once more that in each expression f l5.72f5.73|5.74p all IR diver- 
gencies cancel for arbitrary 6 and only the finite part is left. 

Two comments are in order. First, this decomposition is valid only in the leading 
order in a. In the next orders the inclusive cross-section requires extra emitted particles 
that takes us away from the class of the MHV amplitudes. It is not clear whether in 
this case one has separate IR-safe sets or everything is mixed together and only the total 
cross-section is finite. In the latter case one probably faces the complication since the non- 
MHV amplitudes are not known to possess a simple structure as the MHV ones, though 
the origin of this simplicity is unclear. The second comment concerns the contribution of 
the matter fields. In the leading order we singled it out in the class C. At the same time, 
in general, there is the contribution of them to the virtual part and to the splitting one 
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via the gluon distribution function. They contain the 1/e terms. However, the matter 
field contribution to the virtual part is proportional to the tree-level 2x2 cross-section 
with the coefficient (|?T,y -|- ^Ug) [5Z] and the contribution to the splitting function comes 
with the /3-function, i.e. with the same coefficient but with the opposite sign. Thus, 
these contributions have the same structure and completely cancel each other. So, in the 
leading order our separation becomes possible. 

Defining now the physical condition for the observation we get several infrared-safe 
inclusive cross-sections 



Registration of two fastest gluons of positive helicity 



A 



MHV 



5=1/3 



B 



anti-MHV 



<5=1' 



Registration of one fastest gluon of positive helicity 



MHV 



A 



Anti-MHV cross-section 



<5=l/3 



B 



anti-MHV 



6=1/3 



+ c 



A'latter 



<5=1' 



B 



anti-MHV 



(5=1 



+ c 



Matter 



5=1 



(5.75) 



(5.76) 



(5.77) 



Relative simplicity of the virtual contribution (I4.36P i which contains logarithms and 
no other special functions suggests a i similar structure of the real part. However, this 
is not the case. While the singular terms are simple enough and cancel completely, the 
finite parts are usually cumbersome and contain polylogarithms. The only expression 
where they cancel corresponds to the 5 = 1 case which is possible only for the last set 
of observables, namely, for the anti-MHV cross-section (15.771) . Choosing the factorization 
scale to be Q/ = we get: 



da 



13 / 

a 

477 




(5.78) 



anti-MHV 

2(c^ + 2c3+4c2 + 6c+19)log2 



l-c ^ 
2 ^ 



(1-C)2(1+C)4 
(C2 + I)l0g(i±^)l0g(i^) 



+ 2 



[c*~2c^+4:c'^-6c+19) log^m^) 

(1-C)4(1 + C)2 



67r2(3c2 + 13) - 5(61c2 + 99) 



-2 



(llc^ 



(1- 
-31c2 



C2)2 

-47c- 133) log( 



l-c ^ 
2 ^ 



3(l + c)3(l 



+ 2 



9(1-C2)2 

(llc=^+31c2-47c+133)log( 



l±c\ 
2 ) 



3(1-C)3(1 + C)2 

One can see that even this expression does not repeat the form of the Born amplitude 
and does not have any simple structure. While the dependence on the parameter 
of dimensional reduction is completely canceled, the finite answer, as in the toy model 
example, depends on the factorization scale. This dependence comes from the asymptotic 
states which violate conformal invariance of the Lagrangian. This dependence seems to 
be unavoidable and refiects the act of measurement. Construction of observables which 
do not contain any external scale remains an open question. 
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6 Discussion 



Remarkable factorization properties of the MHV amplitudes accumulated in the BDS 
ansatz (with the so far unknown modification) and duality with the string amplitudes 
via the AdS/CFT correspondence seem to suggest the way to get the exact solution of 
the A/" = 4 SYM theory. However, "to solve the model" might have a different meaning. 
Calculation of divergences and understanding of their structure is very useful but surely 
not enough, it is the finite part that we are really for. The knowledge of the S-matrix 
would be the final goal though the definition of the S-matrix in conformal theory is a 
problem. Even in the absence of the UV divergences there are severe IR problems and 
matrix elements do not exist after removal of regularization. 

The purpose of this paper is to present all the details of the calculation with explicit 
cancellation of the infrared divergencies in properly defined cross-sections in the planar 
limit for A/" = 4 SYM. The main results were summarized in our short letter [58]. We 
do obtain IR safe observables in the weak coupling regime in the next-to-leading order of 
PT which are calculated analytically. The same procedure can also be applied to J\f = 8 
SUGRA P]. 

Unfortunately, our calculation has demonstrated that the simple structure of the am- 
plitudes governed by the cusp anomalous dimension has been totally washed out by com- 
plexity of the real emission matrix elements integrated over the phase space. This means 
that either the N' = 4 SYM theory does not allow such a simple factorizable solution 
or that we considered the unappropriate observables that do not bear the impact of the 
A/" = 4 symmetry. One can obviously see the presence of A/^ = 4 supermultiplet in the 
virtual part but not in the real emission. It would be of great importance to find such 
quantities. 

Another unfortunate feature of inclusive cross-sections is the dependence on the fac- 
torization scale. The experience of QCD, which is very similar to the A/" = 4 SYM theory 
from the point of view of the IR problems, tells us that in inclusive cross-sections the 
IR divergences cancel and one has finite physical observables. However, in QCD one has 
confinement and considers the scattering of the bound states (hadrons, glueballs) rather 
than the individual particles. In this case, one usually factorizes the hard part from the 
soft part introducing the factorization scale. The dependence on this scale is canceled 
between the hard and soft parts contrary to our case where only the hard part is present. 
But in QCD one also has an additional scale. The parton distributions are defined exper- 
imentally at some scale Qo and the dependence on this scale is left. This dependence is 
governed by the same DGLAP equations as the dependence on the factorization scale, so 
from this point of view the situation in QCD is not better than in our case. 

In both the cases one has to introduce some parton distributions which are the func- 
tions of a fraction of momenta and, in higher orders, of momenta transferred. This leads 
to the appearance of an additional scale which breaks the conformal invariance. One 
might think of some observables where this scale dependence is canceled, like the ratio of 
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some cross-sections, etc. We have not found such quantities so far, though the construc- 
tion of such truly conformal observables is of great interest. Probably, they might have 
the desired simple structure. 

There is an interesting duality between the MHV amplitudes and the Wilson loop, 
between the weak and the strong coupling regime [251 EH [29]. Perhaps, it would be 
possible, using the AdS/CFT correspondence, to construct the IR safe observables in the 
strong coupling limit (similarly to what we did here) and to shed some light on the true 
calculable objects in conformal theories. 
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Appendix A. Computation of partial amplitudes 

To calculate the cross-section we need the squared matrix elements summed over helicities 
and color. They can be expressed in terms of the corresponding partial amplitudes [5] 

colors 

= g'-'N:-\Nl-l){^^r E l4'-'"°^Hpi,P.(2),...,P.H)r- (A.I) 

For the massless partial helicity amplitudes it is convenient to use the so-called spinor 
helicity formalism initially introduced in [Hni ED 102] (for a review see [S3])- In this 
formalism the on-shell momenta of every z-th external massless particle = Q is 

represented in terms of a pair of massless commuting spinors A^*^ and A^*'* of positive and 
negative chirality in the following way: 



P?^p2=P;:H^'^)aa = A«A«. (A.2) 

The spinor inner product is defined by: 

e-^^A^A? = (A«A(^-)) = (zj), e^>!{^>!f = [A^A^^')] = [tj], (A.3) 
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thus the complex conjugation of the product is 

mr = m- (a.4) 

The scalar product of the two light-like momenta can be represented in terms of these 
products as 

or equivalently 

= Sij, (A. 6) 

where the standard notation {pi + pj)'^ = Sij is used. 

All the tree-level partial MHV amplitudes can be combined into a single A/" = 4 
supersymmetric expression, first suggested by Nair 

/ n \ 
zAf=4 MHV _ x8 I ^ 

\i=l 



where rjl^'^ are the Grassmannian coordinates, A = 1,...,4 is the S'f/(4)/j fundamental 
representation index. Z:^^'^ mhv jg invariant under SU{4){i transformations of rjl'^'' and 
under the cyclic permutations of momentum labels i. In the product nr=i(^?'^ + 1) 
has to identify i + n with i. The Grassmannian-valued delta function is defined in the 
usual way: 



/ n \ 4-|/n \ / ^ \ Tin 

(T.Kvi''] = n I E A-^l^M E aU^M = 1^ n E i^k) (A.8) 

\j=l / A=l ^ \j=l / \fc=l / ''^ A=li,k=l 

So one can rewrite Z^^^ mhv 

= ^ E m {Im) {ah) {do) (r^fMSf^L^i^i%i%f ) (A.9) 



where 

n 

Vn = X{{l,t + l). (A.IO) 

i=l 

Using the Taylor expansion of Z^^'^ mhv powers of 77^"^^ one gets the sum of ( "'•"^~'^^ )^ 
terms each involving a product of 8 distinct r]\'^\ One can identify the coefficient of the 
product of 8 77's in each term in the expansion with a particular tree component partial 
amplitude. It is very useful to define the following differential operators with the self- 
explanatorily notation: 

= 1, 
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J_^ABC'D _ 



4! driY"^ dri\^^ drjf^ dri\^^ dj]^^^ drif^ drjf^ dvjf''^ ^ 

dm 

1 

q {t)a = -—,^ABCD 



3! dr]\^^ drif'^ drif''^ 



drji 'drj. 



2! dril^^dril^^ 

Taking various combinations of products of these operators one can construct a set of 8-th 
order differential operators. These 8-th order differential operators act as projectors on 
the component partial amplitudes: q~^{i)^ corresponds to the fermion q'^''^ of the Af = 4 
supermultiplet, to qj, A^^(i) to A"^-^, and A^B(i) to A^^- 

For example, the Parke- Taylor n-gluon amplitude can be written as: 

A^^^^\g-g-g^...g^) = r (l)r (2)^+(3)...^+(n)Z;f =^ = {12)'^, (A.12) 

and the squared partial amplitude \A^'^'^'^\g^g^ g^ ...g'^)]'^ then takes the simple form (it 
is implemented that momenta are ordered as Pi,P2,P3, ■■■,Pn) 



\Alr'''\g-g-g^...g^)\' = Alr^^\g-g-g+...g+)Aii^^^\g-g-g^...g^r 

(12)4[12]4 sj 



'12 



(A.13) 



To extract from flA.ip some specific helicity configuration for the MHV amplitude, one 
has to sum over the permutations only in the denominator of flA.13p [55J. So, for example, 
for the Parke- Taylor n-gluon amplitude one has 

\M^^'^^\g~g-g+...g^)\' = g^-'N: s{, ^ ^ • (^-14) 

The anti-MHV amplitudes also needed for our computation can be obtained from 
the corresponding conjugated MHV amplitudes. For example the anti-MHV amplitude 
A^{g~ g~ g^ g~ g^) can be obtained from the MHV amplitude A^{g~^ g~^ g~ g~^ g^) by making 
a complex conjugation. 

Below we present the list of four- and five-point tree amplitudes which are relevant to 
our calculation. The four-point amplitudes are 

A'r\g-g-g^g^) = r (l)r(2)^^(3)^+(4)Zf =^ = (12)^:1, (A.15) 

'4 
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At''\9-9^9-9^) = 9-{l)9H2)g-{3)g^4)Zf=' = (13)^:1, (A.16) 

At'\9-q^9^qA) = g-{l)q^'^{2)g^3)qjmf=' = (12)(14)^:1, (A.17) 

At'\9-qA9U^) = g-{l)qAmHm^'H^)2f=' = (14)(12)=^:^, (A.18) 

^^'^H^^A^'^+Aab) = r(l)A''(2)^+(3)A,B(4)Zf =^ = (A.19) 
For the computation of the real emission we need the five-point tree amphtudes 

I {tree 



^5 



(tree 



{9-9-9^9^9^) = r(l)r(2)^+(3)^+(4)^+(5)Z5^=^ = (12)^:1, (A.20) 

'5 

= (^+(l)^+(2)r (3)^+(4)r (5)^5^=^ ^^^)* = [35]^^, (A.21) 

'5 

(^-^7-^7+^7+^7-) = (^+(l)^+(2)r (3)r (4)^+(5)Zf =^ ^^^)* = [35]^^;^, (A.22) 

' 5 

{9-9-9\^qA) = r (l)^+(2)r(3)^^'+(4)gl(5)Zf =^ ^^^)* = (A.23) 

' 5 

(rr^+g^g^) = (^+(l)^+(2)r(3)gl(4)^^'+(5)Zf =^ ^^^)* = i^^, (A.24) 

'5 

4*''^'^^(p-r^+A''AAB) - r(l)r(2)r(3)A^^(4)A^(5)Z5^=^ ^^^)* = i^ffl^.(A.25) 

'5 

We also provide the list of the Born cross-sections used in Sect. 4. 



^5 



< (tree 
^5 



j^itree] 



A (tree 
^5 



(-9+9) ^,2 /\r2 / „2 _L +2 _L ,,2\ /,^2' 



(A.28) 

(A.29) 
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These cross-sections are written down for the set of momenta {pi,P2,P3,P4) with the 
conservation law P1+P2 = P3 +P4,- In the case of initial splitting, according to (14. 54 p . one 
should use the cross-sections calculated for the set {zpi,p2,P3,Pi) with a new conservation 
law zpi + P2 = P3 + Pi- To get them, one should substitute the modified values for the 
Mandelstam variables f l4.64p into f lA .2611X7311) and multiply the cross-sections by the factor 
+ z — c{l — z)Y which comes from the modified delta function 5^{zpi +P2—P3 ~Pa)- 
The same procedure but with the replacement c ^ — c refers to the pi ^ p2 case. 



Appendix B. Calculation of phase space integrals 

Consider the structure of the matrix elements in detail. First off all it is convenient to 
rewrite the standard three-particle phase space 

in the following form: 



»3 = ^*iPl)j§hsH(P.-ki')j0rTi*(k')j0^A2-mm+m-P,-Pd 

(B.2) 



The integral we are interested in is 

I \M,\''Ssips,k,p,-k)d(f)s, (B.3) 

where the matrix element |A^5p for the five-point amplitude consists of 12 terms with 
identical numerator but different denominators. The typical integrand looks like 

H2 2{{p,+p2fY 



S13S25S35S24S14 {pi - P3Y{p2 - kY{-P3 - ky{P2 - [P4 - k]y{pi - [p4 -k]y 

(B.4) 

Our strategy is to use the on-shell conditions to simplify all the terms in the sum 
so that the integral over d^k can be calculated exactly. For the remaining integrals we 
evaluate the necessary terms of the e-expansion. 

Taking into account the conservation of the momentum Pi + P2 = P3 + Pa and the 
on-shell conditions 

pI = 0, pI = 0, pI = 0, k^ = 0, {p,-k)^ = 

one can rewrite the integrand flB.4p as 

r (Pl,P2)^ 



{Pi,P3){P2,k){p3,k){p2,P4 - k){pi,pi - k)' 



(B.5) 
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where we use the notation (p,k)=pk for the scalar product. 
The next step is to use the partial fraction with respect to k 

J ^ {Pl,P2)^ 1 / 1 ^ 1 ^ 



{PuP3){P3,k){pi,Pi- k) {P2,P4) \{P2,k) (P2,P4-/C), 
(Pl,P2)^ 1 , (Pl,P2)^ 



+ 



(Pl,P3)(P2,P4) (P3, k){p2, k){pi,P4-k) {PuP3){P2,P4) (Ps, k){p2, P4- k){pu p^- k) 

{pi,P2y 1 

(P1,P3)(P2,P4) (P1,P4) - (p4,P4)/2 

f I 1 1 \ 



XP3,k)(p2,k) (P3,k)(pi,p4-k) (P2,k)(pi,p4-k)j 

(Pl,P3)(P2,P4) (Pl,P4) + (P2,P4) - (P4,P4)/2 

f 1 1 1 \ 



iP3:k){p2,P4- k) {P3:k){pi,p4- k) {P2: P4- k){pu P4- k) J 

SO that one gets at most two brackets with momentum k in the denominator. 

In the case of momentum k in the numerator, this procedure also works but with some 
variation. For example, one has 



J 



{Pi,k) 



{P1:P2){PUP3){P3,P4 - k){p2,k){k,p4 - k) 

1 {P3,k)+py2-{p2,k) 

(Pi, P2) (Pi, P3) (P4/2) (P2, k) (P3, P4 - k) 

\ (P3,P4) - (P3,P4 - k)+pl/2 - {p2,k) 

(Pl,P2)(pi,P3)(P4/2) {P2,k){p3,P4 - k) 

(P3,P4) +P4/2 

(P1,P2)(P1,P3)(P4/2)(P2, k){p3,P4 - k) 

1 1 



(B.7) 



(P1,P2)(P1,P3)(P4/2)(P2,A;) (P1,P2)(P1,P3)(P4/2)(P3,P4 " k)' 



Since we usually have (p^, kY in the numerator this procedure has to be applied several 
times. This way we increase the number of terms in the integrand but drastically simplify 
the integration. 

The resulting integrals over k have the standard form 

j d^M+((p4-A;)2)5+(A;2)T, , (B.8) 

where 
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Ti 



T, 



1 

1 

(Pi,P4 - kY{pj,pi - kf 



and can be calculated by the method of unitarity. They correspond to the box-type 
diagrams and one can perform the cuts and then take the imaginary part. For example, 
the integral 



{pi + kY{p2 - pa + ky 



(B.9) 



where Pi= P2= pi = ^ and P4 7^ can be obtained from the box-diagram shown in FigJHl 




^4 ^4 



Figure 6: The box diagram corresponding to the integral fIB.Qp 



For the first time this integral was calculated by van Neerven [65] and the answer is 
given by 



(p1/4)- 



0(P4) 



X 27r 2 



(PiPiYiPjPAf 

_,_,T{D - 3)r(D/2 - 1 - anD/2 - 1 - b) 



(B.IO) 



D 



T^{D/2-l)TiD-2-a-b) ^^^(^^^'^2 



where 



Ti 
T2 
T3 



1 - 



'^{Pt,P4){Pj,PA)' 
{Pi,Pj){P4,P4) 
2(Pi,P4)(Pi,P4)' 
^ _ {Pi,Pj){P4,P4) 
2(Pi,P4)(Pi,P4)' 
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Removing the integral over d with the help of the delta function we are left with 
the last integration over d^p^. Using 5'^ (pi) one can take the integral over pg and going 
to the spherical coordinates d^~^p^ — |p3|^~^(i|p3|dfii3 arrive to the single integration 
over the modulus of |p3|. 

Here we face the problem of singularity at IP3I =0. It comes from the delta function 
in the integration over pg and in some cases is not compensated by the matrix element. 
For two matrix elements corresponding to the MHV {g'^g'^ — > g'^g'^g'^) and anti-MHV 
{9'^ 9^ ~^ g^g^g~) amplitudes, the first case is singular while the second is not. However, 
as we explained earlier, we cut the integral over \p^\ at (1 — S)E/2 and no singularity 
appears. 

Let us now turn to the calculation of the last integral. Since ps is a dimensionful 
parameter, it is appropriate to change variable to dimensionless one using 

|P3l = f(l-^)- (B.ll) 

Then the integral over x goes from to 6. 

The typical integral to be calculated is of the form 

/ dx ^ ^ 2Fi{l,-e;l-e;qx^il-xr), (B.12) 

Jo 1+ px 

where m and n take the values 



(m,n) = {(l,0), (0,2), (1,-2), (-1,2)}. 

For our purposes we need to calculate this integral to the order 0{e). The source of 
divergence is the singularity at a: = 0. When 6^1, one can expand the hypergeo metric 
function in e up to the order e and then calculate the integral. Then for the configurations 
(1,0), (0, 2) and (1,-2) the calculation is straightforward while for the case of (—1, 2) one 
first makes the transformation of the argument of the hypergeometric function from z to 
1/z 

2Fi{a,b;c\z) = r(c - a)r{b) ^ ' [a,a - c+ l;a - b + (-^)-" 

and then apply the expansion. 

For example, consider the integral 



/ 

Jo 



s 

dxx~^-'{l - x)-^'2Fi{l, -e; 1 - e; qx-'^{l - x)^). (B.14) 
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After applying the transformation flB.lSp it is reduced to the following form 



dx\T(l-e)T(l + e) 



' q{l — x) 



X 



+ 



xT{-l - e)r(l - e)2F^ (l, 1 + e; 2 + e; 



g(l -x)2r(-e)2 
Performing the expansion over e one gets to the order of 0{e) 



(B.15) 



l2 = - 

-log 



+ 



log 5 + log q log (5 — log { (5 — 1) q) log 5 



2(,5-l)g- y/T^+l 
2g 



+ log (q5^ -2q5 + S + q) log S 



log (5 — log 



2((5-l)g + v/r^4^+l 



log g or ■ n x\ ^ 
— 2Lt2{l-6) 



2q 

2 



log 5 



12 



-log 
-log 
+ log 
+ log 
+Li2 



-2g+ v/r^4^+ 1 
2q 

2q \ 



j log 
log 1^ 
log 



2q 



-2g + v/T^4g+l 



^2g + ^1^4^-1, 
2^g 

^2g + v/r^-l^ 
2(5g 

"-2g + ^r^+l 
'2{6-l)q + ynr^+l 
-2g + v/T^4^+l 



2g + yr^-n 

2g J 
2(5-l)g-yr^+l~ 



2g 



log 



2(5-l)g + yr^4^+l 



2g 



Li'2 



-25q + 2q + ^TT^ - 1 
2q + v/r^4^ - 1 



e. 



Note the singularity when 5 — ?■ 1 in this expression. 

The case of 5 = 1 is more tricky. Here one has the overlapping of two singularities. 
The argument of the hypergeometric function goes to the edge of the circle of convergence 
and it is convenient to use the integral representation 



1 C — JO 



r(c - b) 

As a result, one has a two- fold integral 



ffl^^f,^ /' dx dt t''-\l - ty-'-^x'^-Ul - x f'\l - tqx"\\ - 
r(c — b) Jo 



(B.16) 



(B.17) 



where parameters a, b and c take the values a = 1, b = — e, c = 1 — e. Choosing particular 
values of a,(3,m,n one can observe the overlapping divergencies. Consider, for example. 
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the integral 

Cdxdtr^-'x-^-Hl-x)-^^ = Cdxdtt-^-'x-\l-x)-'^' }- (B.18) 

Jo ^ ' l-qtik^ Jo ^ ' x-qtil-xY' ^ ' 

where we see in the last term that the denominator equals zero at t = and x = 0. The 
divergence, which occurs in this case, is the overlapping IR divergence and to handle it 
we use the following trick: we insert in the integral a unity 

1 = 0(a;-t) + 0(t-x), 

which splits the integral into two parts. The first 6'-function gives 

/ dxdz z-^-'x-^'^'{l - ]- -r, (B.19) 

1 — qz[l — xY 



while the other leads to 

r-l 







dz dt r^-~^'z~'(l - zt)-^' -. (B.20) 

^ ^ z-q{l-zty ^ ' 



The calculation now is straightforward. One has to extract a few terms of the e 
expansion. 

For example, the first three terms of the e-expansion for the integral f lB.14p when 
5 = 1 are 
"1 

dxx-^-'{l - xy^'iFiil, -e; 1 - e; qx'^^l - xf) (B.21) 
_i _ ( I ^ ] , _ f- _^ ^ ] .^^1, 0(.=). 



Appendix C. Splitting functions 

The splitting functions Pij which we use to calculate the splitting contribution to the 
cross-section can be obtained from the coUinear limit of the color ordered tree level par- 
tial amplitudes. Suppose one has an n-point partial tree amplitude in < A/" < 4 
supersymmetric gauge theory 

\Pa{l)^ ■■■1 Pa{i)i ■■■1 Pa{n))i 

where a{i) is the color index of i-th particle and Aj is it's helicity. 

It can be shown [5l [6] that the MHV amplitudes have the following universal behaviour 
in the coUinear limit when momenta of two particles i and i + 1 become collinear + 1 

Y.Spl^t-x{a'\b''^\z)A':ll{...,p',,...), (C.l) 

A.c 
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where the two momenta satisfy 

Pi = zp, pi+i = (1 - z)p, 

p being some arbitrary momentum. The sum goes over all possible helicities and particle 
types for which A*^^^ is nonvanishing. The function Split_x{a^\b^^+^ , z) depends on p 
and z. Notice the flip of helicity in Split _x{(^^\ h^^+^ ■, z) which comes from considering all 
particles as outgoing ones. 

Then the polarized version of the splitting function Pij can be obtained from Split 
up to the terms proportional to 6{1 — z) by means of 



= {p^+Pi+lf\Svl^t-x{a'%h'^^\z)\' (C.2) 

and corresponds to the process c — )■ + 1 when the particle c with momentum p and 
helicity A splits into coUinear particles i and i + 1 with momenta zp and (1 — z)p and 
helicities A, and Aj+i, respectively. 

For example, the splitting function Pg+g+ can be obtained from the partial amplitude 

j[iree^g-g-g+g+g+^ ^^^J^^^g ^^^^ Hmit 4||5 = Zp, ^5 = (l " Z)^) ill 



One has 



1 1 (12)4 



Thus, the only one term in the sum (IC.ip survives and A'^^\ in this case is A^l'^'^{g g g^g^). 
This gives 

SpUt.{g\g\z) = ^ , / , (C.4) 
(45) ^z{l - z) 

so that, according to flC.2p . 

All the splitting functions necessary for our computation can be obtained in a similar 
fashion. They look like 



p9 

^9+9- 



Z (l-z) 

3 



Z 



1 
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Pl^ = z{l-z) 



The contributions proportional to 5(1 — 2;) are calculated separately from the requirement 
of conservation of momenta and are absent in our case since they are proportional to the 
/3 function which vanishes in the J\f = 4 SYM theory. 

The "plus" prescription in the expression ^^^J-^^^^ in (1C.6P should be understood in the 
usual way: 



dz- — = / dz—- — . (C.7) 

(1 — zj+ Jo (1 — z) 

When f{z) contains the theta function like in the splitting counterterm 

f{z) = Q{Z - Zmin)9iz), 

one has 

' dZyM_ = r - ZrnMz) - 9(1) 

(1 - z)+ Jo (1 - z) 



,1- Z) Jo 1 - Z Jz^.r^ (1 - 2) 

The splitting function Pqq{z) fl3.23p used in our toy model example can be obtained 
from the polarized splitting functions 

1 



p: 



1 ' 



(1-^)+' 

by summation over helicities. The term proportional to S{1 — z) is obtained from the 
requirement of conservation of the number of quarks 



1^ dz q{z, Q}/fi^) = 1 ^ dzPggiz) = 0. 



Appendix D. Finite parts of amplitudes. 

In general all finite parts have the following structure: 

I'inite part = _^ ifsym{c, 6) + {fAsym{c, 6) + fAsym{-C, 6))] , 

where the functions fsym and fAsym contain Lo^ and Polylog functions of c and 6. Below 
we present the expressions for fsym{c,6) and fAsym{c,S)- 

^^To make the expressions more compact we use L for the logarithms 
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fin 



fsym {C:0) = (1^^5)2 4(13 + 3c )L(5)L(1 - 5) (D.l) 



+ 10(3+c2)L2(<5) - 4(5 + c2)L(^)L(^) - 16(3 + c2)L2(1-<5) - A{9c^ + 35)Lt2{5), 
where 

5i = ^(3 + c2)7r'(l-(5)2 + 32(4-35)5, 

^2 = 4(3c^(l- 5)^ + 37- 265 + 50, 
^3 = -45(0^(5-1) + 115 -15). 



- 2(-l + 4c + c')L2(i^) - 8(3 + 2c + c2)L(l - S)L{^^) 

+ 16cL(5)L(l^)+4(l + c)^L(l^)L(i±A^|li^) 

+ 4(1 + c)^L(l - 5)L( ^ + ^-f-^) )-4(l + c)-L(5)L( ^+^-f-^) ) 

+ 8(l + c)L^2(^)+4(5 + 2c + c2)L^2( ~ ^ 



1 + c ^ 

4(5 + 2c + c')Li2(^^ — ^'-^ — ^), (D.2) 



where 



A = 4(c2(l-5)2-6c(l-5)2 + 5 + 25-352), 

^2 = -4(c2(l -5)2 -2c(3-45 + 52) + 5 + 25-35=^) 



, \(— ++-)\ 

, general b. 

d^^^ J Real J f,^' ^ 



fsym'^'\c,S) = S,+S2m+SsL{l-S) + 10{c' + 3)L\6) (D.3) 

o^Q7^i8 2^ 4^,. , 8(6 + 9c2 + c^) ^^l-c^^^l + c^ 
- 2(37 + 18c +c )Lz2(5) + (iZT^^ 
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where 

8{3c^ + 5)6^ + 3{7c^ -95)6^ + 6{67c^ + 513)6 2.^ 64(llc2 + 7) 

^ 9 + 3^' + 3(c2-l) 

32(c3 + 12c2 + 19c + 9) 32(2c=^ - 120^ + 19c - 9) 



+ 



3(1 - c)(l + 5 + c(l - 5)) 3(1 + c)(l + 5 - c(l - 5)) 
32(c3 + 4c2 + 5c + 2) 32(c3 - 4c2 + 5c - 2) 



+ 



2' 



3(l-c)(l + 5 + c(l-5))2 3(l + c)(l + 5-c(l-5)) 

3 c"^ — 1 3(c^ — Ij 

32(c3 + 5c2 + 11c + 7) 32(c3 - 5c^ + 11c - 7) 

^ (1 - c)(l + 5 + c(l - 6)) " (1 + c)(l + 5 - c(l - S)) 

32(c3 + 4c2 + 5c + 2) 32(c3-4c2 + 5c-2) 



+ 



(l-c)(l+5 + c(l-5))2 (l + c)(l+5-c(l-5))2 

64(c''^ + 3c2 + 3c + 1) 64(c=' - 3c^ + 3c - 1) 



3(1 - c)(l + S + c{l- <5))3 3(1 + c)(l + 5-c{l- 5)) 



3' 



8/0 2 .^f3 /4 o2 ^.^:2 32(c^+8c2- 11)5 3c6+251c^+2953c2+313 

--(3c^ + 5)53_(,4^2,2_75),2_^___^^ 

61(c-' - 5c- + 11c - 7) 64(c^ + 5c- + 11c + 7) 



+ 
+ 



(l + c)(l + (5-c(l - 




(l-c)(l + (5 + 


c(l- 


-s)) 


64(c3 - 4c2 + 5c - 


2) 


64 (c-*^ + 4c2 H 


h5c 


+ 2) 


(l + c)(l + 5-c(l- 


■ s)r 


(l-c)(l + 54 


cil 




128(1 + c)3 




128(1 

+ rr—— 


-c) 


3 



3(l-c)(l + 5 + c(l-5))3 3(l + c)(l + 5-c(l-5)) 



^ 1|±|£.(1^).(1±AZ|(1^) p.4) 

4(7 + c^)(l-c) l + c l + 5-c(l-5) 

(TT^) ^^^^^^ 2 ^ 

8(c4-12c3 + 34c2 + 12c-43)^,^ A + 5 - c{l - 5), 

OTcj-' - ^^^^ 2 ) 

2(c^-2c3 + 8c^-6c+15) ^2^ 1-^ ^ 

. r/l -C^ , ^ ,l + 5-c(l -5), 
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4 (3c^ - 12c^ + 46c2 + 12c - 33) 1 + 5 - c(l - 5) 

{iTcP ^^^^^^ 2 ^ 

(c^ - + 3c^ - 76c=* - 153c2 + 14c + 149) ^ , 

+ iT^, -i^'^i-—/^ 



where 



A = 7 T (2597 + 2405 - 105^2 + 245^ + 3(:^(6^ - 1) + 3c^(3 + 6^) 

6(1 + c) V \ J \ J 

-2c3(lll - 24(5 + 95^ - 46^) + 2c2 (489 +1205 -695^+20(52) 



-c(2655- 2405+2255^ -565^)) , 



A2 = , (3(5' - + 3(5' + 3)c^ + 2(85=^ - 35' + 365 - lll)c2 



6(l + c 

+6(85^ - 175' + 285 + 163)c' + 3(165^ - 635' + 1045 - 885)c 
+165^ - 935' + 2165 + 2597) . 

In the case 5 = 1 one gets major simphfications: 



(__++_), _ 22 57 - 937r'-3c^7r' + c' (303 - 487r') 

JSym (C, Ij — - 



where 



Ai = -^(3c^-46c2 + 280c' - 646c + 689), 

A2 = -^(c^ + 2c^ + 7c^ + 68c2- 121c' -38c + 209), 
Az = c^ - c^ - 6c2 - 178c' - 603c - 493. 



1^- 32(c' + l)% l-c l + c 47r'(l-c-^) + 132(c' + 3) 
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where 



Ai = -(3c^-44c^ + 222c2- 450c + 277), 

A2 = -2(c^ + 2c^-2c^ + 50c-67)(l-c), 
^3 = -4(c^-2c=^-18c2-146c-211)(l + c). 



4r^'^(c, 1) = -24(c^ + l)L(l^)L(l±^) + 6(llc^-3) - 7^^(l-c^)^ (D.9) 



2 (3c^ - 47c=^ + 213c2 - 369c + 184) ^ . 1 - 

TT^ ^^^^ 

3 (c^ - 3c2 - 17c - 125) (1 + c)2 1-c, 
+ ^ (l-c) ^ ^M— ). 



d(7- 



2^3 



tin 



general S. 



^^^13 JinSplit J fi 

fty;n'''"\c,S)^S,L{l-S)+S2, (D.13) 



44 



where 

16,2 .^.3 o/ 2 r^.2 ../2 . ^^ . 8(270^ + STSc^ + 59) 

5i = — {0^ + 1)6' + 8{c^- 5)6^ + 16{c^ + 17)6-^ 3(^^) 

64(c=^ + 5c2 + llc + 7) ^ 64(c=^ - 5c^ + 11c - 7) ^ 64(c=^ + + 5c + 2) 



(l-c)(l + 5 + c(l-5)) (l + c)(l + <5-c(l-5)) (l-c)(l + 5 + c(l-5)2 
64(c3 - 4c2 + 5c - 2) 128(c''' + 3c2 + 3c+l) ^ 128(c''^ - Sc^ + 3c - 1) 



(l + c)(l + 5-c(l-(5))2 3(l-c)(l + (5+c(l-5))3 3(l + c)(l + 5-c(l-5))3' 

16,, xc-^ 4, „ 8, . 256(2c2 + l) 

52 = --(c^ + 1)5^ - -(5c2- 13)5^ --(llc^ + 89)5-^^^^-^ 

32(2c3 + Qc^ + 12c + 5) 32(2c3 - Qc^ + 12c - 5) 
~ 3(1 - c)(l + 5 + c{l- 6)) ^ 3(1 + c)(l + 6- c(l - 6)) 
32{c^ + 3c2 + 3c + 1) 32(c3 - 3c^ + 3c - 1) 

^3(l-c)(l + 5 + c(l-5))2 ~ 3(l + c)(l + 5-c(l-5))2- ^ ^ ^ 

(__++_) _ 16(1 - c)(4c2 - 17c + 37) / l + ^-c(l-^) ,4-c\ 

+ 4(3.c^)(i^(^^:^)-^(V))^ 

+ ^^-^-^^f+f--^^)L(l-^)L(^ + ^-f-^)) (D.15) 
+ 8{3 + c'')(li2{-I^S)-Li,' '^^ 



1 + c ' 'l + 5-c(l-(5)7 

8(c=^ + 15c2 + 51c + 45) f^. ^1-c^ 1 \ 

+ - Lt2{-{1 - c){l - 6)) j . 

In the case 5 = 1 one gets major simphfications: 

1) = ^ ((67r^ - 49)c2 + 187r^ - 415) , (D.16) 

16(l-c)(4c2-17c+37) ^^l-c^ , 16(l+c) (c2 + 6c+21) ^ . ^1-c^ 

= 3(rT^^ ^[-^) + 1^^ H2 



V «"i3 JinSvlit / 



(D.17) 

16, 

y 

64(4c^-17c+19)(l-c) ^,l-c, 64(c + 3)^(l + c) 1 - c 

fAsym {C, 1) - ^^^^-^ L(^) ^^2(^). 

(D.19) 

45 



fk^''''\c,l)^4i^c' + 23), (D.18) 



^\ '^^is ) InSplit 




/g„;™(c,l) = -16(3c^-l), (D.20) 

16(13 - 4c)(l - c)%^l -c^ ^ 48(c + 5)(l + c)% . ^1-c^ 
1 + c) 



f Asym (C, 1) = 7^^:^ H^^) + n ' Lt2{^^) . (D.21 
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